Inverted Berezinskii-Kosterlitz-Thouless Singularity and High- Temperature Algebraic 
Order in an Ising Model on a Scale-Free Hierarchical-Lattice Small- World Network 
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We have obtained exact results for the Ising model on a hierarchical lattice incorporating three key 
features characterizing many real-world networks — a scale- free degree distribution, a high clustering 
coefficient, and the small- world effect. By varying the probability p of long-range bonds, the entire 
spectrum from an unclustered, non-small-world network to a highly-clustered, small-world system 
is studied. Using the self-similar structure of the network, we obtain analytic expressions for the 
degree distribution P(k) and clustering coefficient C for all p, as well as the average path length 
£ for p = and 1. The ferromagnetic Ising model on this network is studied through an exact 
renormalization-group transformation of the quenched bond probability distribution, using up to 
562,500 renormalized probability bins to represent the distribution. For p < 0.494, we find power- 
law critical behavior of the magnetization and susceptibility, with critical exponents continuously 
varying with p, and exponential decay of correlations away from T c . For p > 0.494, in fact where 
the network exhibits small-world character, the critical behavior radically changes: We find a highly 
unusual phase transition, namely an inverted Berezinskii-Kosterlitz-Thouless singularity, between 
a low-temperature phase with non-zero magnetization and finite correlation length and a high- 
temperature phase with zero magnetization and infinite correlation length, with power-law decay 
of correlations throughout the phase. Approaching T c from below, the magnetization and the 
susceptibility respectively exhibit the singularities of exp(— C/y/T c — T) and exp(D/^/T c — T), with 
C and D positive constants. With long-range bond strengths decaying with distance, we see a phase 
transition with power-law critical singularities for all p, and evaluate an unusually narrow critical 
region and important corrections to power-law behavior that depend on the exponent characterizing 
the decay of long-range interactions. 

PACS numbers: 89.75.Hc, 64.60. Ak, 75.10.Nr, 05.45. Df 
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I. INTRODUCTION 



Complex networks provide an intriguing avenue for 
tackling one of the long-standing questions in statisti- 
cal physics: how the collective behavior of interacting 
objects is influenced by the topology of those interac- 
tions. Inspired by the diversity of network structures 
found in nature, researchers in recent years have investi- 
gated a variety of statistical models on networks with 
real- world characteristics [H S 01 Three empirically 
common network types have been the focus of atten- 
tion: networks with large clustering coefficients, where 
all neighbors of a node are likely to be neighbors of each 
other; networks with "small-world" behavior in the av- 
erage shortest-path length, £ ~ ln(iV), where N is the 
number of nodes; and those with a power-law (scale- free) 
distribution of degrees. Since the pioneering network 
models of Watts-Strogatz Q , which exemplified the first 
two properties, and Barabasi- Albert p|, which showed 
how the third could arise from particular mechanisms of 
network growth, significant advances have taken place in 
understanding how these properties affect statistical sys- 
tems. The Ising model has been studied on small-world 
networks H 0,11, HG3, along with the XY model [HI, 
and on Barabasi- Albert scale- free networks 0, 0] . On 
random graphs with arbitrary degree distributions, the 
Ising model shows a range of possible critical behaviors 
depending on the moments of the distribution (or in the 



specific case of scale-free distributions, the exponent de- 
scribing the power-law tail) |l4l llfil| , a fact which is ac- 
counted for by a phenomenological theory of critical phe- 
nomena on these types of networks |l6j |. 

In the current work we introduce a novel network struc- 
ture based on a hierarchical lattice 0,0,0 augmented 
by long-range bonds. By changing the probability p of 
the long-range bonds, we observe an entire spectrum 
of network properties, from an unclustered network for 
p = with £ ~ N 1 / 2 , to a highly-clustered small- world 
network for p — 1 with £ ~ IniV. In addition, the net- 
work has a scale-free degree distribution for all p. Due 
to the hierarchical construction of the network, together 
with the stochastic element introduced through the at- 
tachment of the long-range bonds, this network combines 
features of deterministic and random scale-free growing 
networks ^ El El El |U |H |H ^ , and in the 
p = 1 limit its geometrical properties are similar to the 
pseudofractal graph studied in Ref . [2^ • The self-similar 
structure of the network allows us to calculate analytic 
expressions for the degree distribution and clustering co- 
efficient for all p, as well as the average shortest-path 
length £ in the limiting cases p — and 1 . 

A renormalization-group transformation is formulated 
for the Ising model on the network, yielding a variety 
critical behaviors of thermodynamic densities and re- 
sponse functions. For the quenched disordered system 
at intermediate p, we study the Ising model through 
an exact renormalization-group transformation of the 
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quenched bond probability distribution, implemented nu- 
merically using up to 562,500 renormalized probabil- 
ity bins to represent the distribution. We find a finite 
critical temperature at all p, with two distinct regimes 
for the critical behavior. When p < 0.494, the mag- 
netization and susceptibility show power-law scaling, 
and away from T c correlations decay exponentially, as 
in a typical second-order phase transition. The mag- 
nitudes of the critical exponents, which continuously 
vary with p, become infinite as p — > 0.494 from be- 
low. For p > 0.494, in fact coinciding with the on- 
set of the small-world behavior of the underlying net- 
work, we find a highly unusual infinite-order phase tran- 
sition: an inverted Berezinskii-Kosterlitz-Thouless sin- 
gularity HIES, 

between a low-temperature phase with 
non-zero magnetization and finite correlation length, and 
a high-temperature phase with zero magnetization and 
infinite correlation length, exhibiting power-law decay 
of correlations (in contrast to the typical Berezinskii- 
Kosterlitz-Thouless phase transition, where the alge- 
braic order is in the low-temperature phase). Approach- 
ing T c from below, the magnetization and the suscep- 
tibility respectively behave as exp(— C/\fTc — T) and 
exp(D/\/T c — T), with C and D calculated positive con- 
stants. 

Infinite-order phase transitions have been observed for 
the Ising model on random graphs with degree distribu- 
tions P(k) that have a diverging second moment (k 2 ) [Til 
Hfij , but for these systems T c — oo on an infinite network. 
An infinite-order percolation transition has been seen in 

models of growing networks HJIHIHHIslSIIllii, 
with exponential scaling in the size of the giant compo- 
nent above the percolation threshold. A prior observation 
of a finite-temperature, inverted Berezinskii-Kosterlitz- 
Thouless singularity similar to the one described above 
was made in Ising models on an inhomogeneous growing 
network |39| and on a one-dimensional inhomogeneous 
lattice [MM El El- 

The final aspect of our network we investigated was 
the effect of adding distance-dependence to the interac- 
tion str eng ths of the long-range bonds, along the lines 
of Ref . [l(j , where distance-dependent interactions were 
considered in a small-world Ising system. With decaying 
interactions, the second-order phase transition for all p 
has a strongly curtailed critical region and corrections to 
power-law behavior that vary with the exponent a de- 
scribing the decay of interactions. 

This paper is organized as follows: In Sec. II, we de- 
scribe the basic geometric features of the network, start- 
ing with the construction of the lattice (II. A), whose 
deterministic nature allows us to derive exact results 
on the degree distribution (II.B.l), clustering coefficient 
(II. B. 2), and average shortest-path length (II. B. 3). Ad- 
ditional details about the derivations in this section are 
included in Appendix A. In Sec. Ill, we study the Ising 
model on the network through exact renormalization- 
group transformations, first for the simpler cases where 
the long-distance bonds are either all absent (III. A), or 



all present (III.B). In the latter case, we look at both 
uniform interaction strengths along the long-distance 
bonds (III.B. 1) and interactions decaying with distance 
(III.B. 2). Finally, we turn to the more complex situation 
where there is a quenched random distribution of the 
long-distance bonds in the lattice, which requires formu- 
lating a renormalization-group transformation in terms 
of quenched probability distributions (III.C.l). Analy- 
sis of the flows of these distributions gives a complete 
picture of the thermodynamics of our system over the 
entire range of p (III.C.2). We present our conclusions in 
Sec. IV. 



II. HIERARCHICAL-LATTICE SMALL- WORLD 
NETWORK 

A. Construction of the Lattice 

We construct a hierarchical lattice [l7llT8LIT9l | as shown 
in Fig. The lattice has two types of bonds: nearest- 
neighbor bonds (depicted as solid lines) and long-range 
bonds (depicted as dashed lines). In each step of the 
construction, every nearest-neighbor bond is replaced ei- 
ther by the connected cluster of bonds on the top right of 
Fig. ^with probability p, or by the connected cluster on 
the bottom right with probability 1 — p. This procedure 
is repeated n times, with the infinite lattice obtained in 
the limit n — > oo. The initial (n = 0) lattice is two sites 
connected by a single nearest-neighbor bond. An exam- 
ple of the lattice at n — 4 for an arbitrary p ^ 0, 1 is 
shown in Fig. [21 

The p = case, with no long-range bond, is the hierar- 
chical lattice [l^l on which the Migdal-Kadanoff [2j| 
recursion relations with dimension d = 2 and length 
rescaling factor 6 = 2 are exact. As will be seen be- 
low, the network in this case exhibits no small-world fea- 
ture, with a clustering coefficient C = and an average 
shortest-path length I that scales like A 1 / 2 , where N is 
the number of sites in the lattice. The p — 1 case, on the 
other hand, shows typical small-world properties, with 
the presence of long-range bonds giving the high clus- 
tering coefficient C = 0.820 and an average path length 
which scales more slowly with system size, £ ~ In N. By 
varying the parameter p from to 1, we continuously 
move between the two limits. These and other network 
characteristics of our hierarchical lattice are discussed in 
detail in the next section. 



B. Network Characteristics 

1. Degree Distribution 

After the nth step of the construction, there are a total 
of N n = | (2+4™) sites in the lattice. We categorize these 
sites by the number of nearest-neighbor bonds attached 
to the site, fc nn , and the maximum possible number of 



3 




with probability p 



with probability 1 -p 



FIG. 1: Construction of the hierarchical lattice. The solid 
lines correspond to nearest-neighbor bonds, while the dashed 
lines are long-range bonds, which occur with probability p. 




FIG. 2: An example of the hierarchical lattice after n = 4 
steps in the construction, for p — 0.6. 

long-range bonds attached to the site, k\d, of which on 
average only pk\& actually exist. At the rath level there 
are 4"-™ +1 / 2 sites with fc nn = 2 m , fc ld = 2 m - 2, for m = 
1, . . . , n. In addition, there are two sites with fc nn = 2", 
fcid = 2™ — 1. Thus, the non-zero probabilities that a 
randomly chosen site has degree k are 



Pn{k) = 




respectively for 

k =2 m + r, < r < 2 m - 2, 1 < m < n - 1, 
k =2™ + r, < r < 2" - 2, 

k =2" +1 - 1. (2) 

Since the degree distribution is not continuous, the 
exponent describing the power-law decay of degrees is 
extracted from the cumulative distribution function 
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in the n —*■ oo limit, P C um(fc) = 2fc'=fc -^(^Oj w here 
P(k) = linVj^oo P n (k). For a scale- free network of ex- 
ponent a, Pcum(k) ~ fc 1_Q . In our case P cum (k) ~ fc~ 2 
for large /c, giving a — 3, a value comparable to the ex- 
ponents of many real- world scale- free networks [l] . The 
maximum degree fc mM in the scale-free network should 
scale as fc max ~ Nn a ^ 01 1 which is indeed satisfied, 
for large n, in our network. The average degree (k) n after 
n construction steps is 

{k) n = £ kP n (k) = 3 + p - , (3) 

fe=l 

which goes to (fc) = 3 + p in the infinite lattice limit. 

2. Clustering Coefficient Cm 

If a given site in the network is connected to k sites, 
defined as the neighbors of the given site, the ratio be- 
tween the number of bonds among the neighbors and 
the maximum possible number of such bonds k(k — l)/2 
is the clustering coefficient of the given site The 
clustering coefficient C of the network is the average of 
this coefficient over all the sites, and can take on values 
between and 1, the latter corresponding to a maxi- 
mally clustered network where all neighbors of a site are 
also neighbors of each other. For our network in the 
ii — > oo limit, C can be evaluated exactly: The frac- 
tion limn^oo A n - m+1 / 2 /N n = 3 • 4T m of the sites, with 
k nn = 2 m and k^ = 2 m — 2, have the average clustering 
coefficient C m , where C\ = p and C m for m > 1 is, as 
derived in Appendix A.l, 

r=0 r'=0 \ / \ / 

2^'(1 -„)2 m -2-r-r' | 2r + pfY) f^} 

(2 m + r + r') (2 m + r + r'-l) ' ^ 

We plot the clustering coefficient C 

oo 

C = Y, 3 • i- m C m . (5) 

m— 1 

as a function of p in Fig. [3] Note that C increases almost 
linearly from at p = to 0.820 at p = 1, as can also be 
seen from the expansion of Eq. JSJ to second order in p, 

C = 0.837p- 0.0378p 2 + O(p 3 ). (6) 



3. Average Shortest-Path Length £ n 

Let dij be the shortest-path length between two sites i 
and j in the network, measured in terms of the number of 
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Stage in lattice construction n 



FIG. 3: Clustering coefficient C for the infinite lattice as a 
function of the probability of long-range bonds p. 



bonds along the path. The average shortest-path length 
£ n is the average of dij over all pairs of sites i,j at the 
nth level. For general p we have evaluated this quantity 
numerically. For p — and p = 1 we have obtained exact 
analytic expressions (Appendix A. 2), revealing qualita- 
tively distinct behaviors: For p — we find 

2 n (98 + 27 • 2" + 42 • 4" + 22 • 16" + 21n ■ 4") 



21(2 + 5 • 4" + 2 • 16 r ' 



11 
21' 



2™. 



(7) 



1 /2 

and since N n ~ 4™ for large n, we have £„ ~ N n ' . Com- 
paring this result to that of a hypercubic lattice of di- 
mension d, where the average shortest-path length scales 
as we see that £„ for the p = network has 

the power-law scaling behavior of the square lattice. For 
p = 1, on the other hand, we find 



23 + 4 ■ (-2)" + 44 ■ 4" + 10 ■ 16™ + 6n ■ 4" + 12n • 16" 



9(2 + 5 • 4" + 2 • 16") 



2n/3 , 



(8) 



which means that £ n ~ ln(7V„) for large n. This much 
slower, logarithmic scaling of £ n with lattice size, together 
with the high clustering coefficient, are the defining fea- 
tures of a small- world network. 

In Fig. 0] we show £ n calculated for for the full range 
of p between and 1, for n up to 6. It is evident that 
even a small percentage of long-range bonds drastically 
reduces the average shortest-path length, and that l n 
shows small-world characteristics, scaling nearly linearly 
with n, for p > 0.5. We shall see below that the small- 
world structure at larger p translates into a distinctive 
critical behavior for the Ising model on this network. 



FIG. 4: Average shortest-path length £ n for level n, shown 
for various values of p between and 1. For p = and p = 1, 
£ n is given exactly by Eqs. Q and @. For other p, we have 
calculated £ n numerically, with an accuracy of ±0.3%. 



III. ISING MODEL ON THE NETWORK 

We study the Ising model on the network introduced 
in the previous section, with Hamiltonian 

—pH =J SiSj + KijStSj 

(ij')nn (ij')ld 

+ H B E ^ + S J') + Hs H Sl ' 

<y)„„ » 

where J, > 0, (ij) n n denotes summation over nearest- 
neighbor bonds, and (zj)id denotes summation over long- 
range bonds. We generalize the above, by introducing a 
distance dependence in the interaction constants K%j , 



Kij = Jm y ° 



(10) 



Here the exponent a > 0, and my measures the range 
of the long-range bond between sites i and j: For a lat- 
tice constructed in n steps, those long-range bonds that 
appear at the nth step have my = 1, those that appear 
at the (n — l)th step have my = 2, and so on until the 
long-range bond that appears at the first step, which has 
rriij — n. The long-range term in the Hamiltonian can 
be rewritten as 



KijS iSj A i X! 



-Ko 



E 



(ii) 



<y>i 



where K q = 



Jq~ a and (ij)id,q denotes summation over 
long-range bonds with my = q. 

The Hamiltonian of Eq. © includes two types of mag- 
netic field terms, one counted with bonds (Hg) and the 
other counted with sites (Hs). We shall calculate the as- 
sociated spontaneous magnetizations at Hb = Hs = 0, 



Mr 



— E 

^ V nn 



(12) 
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where -/V nn = 4" is the number of nearest-neighbor bonds 
after the nth construction stage, so that N nn /N n = 3/2 in 
the limit n — > oo. For a translationally invariant lattice, 
where each site has the same degree, Mb and Ms would 
be simply related by Mb — 2Mg, but for the hierarchical 
lattice this is no longer true due to the different degrees 
of the sites. 

Before turning to the phase diagram and critical prop- 
erties of the system for general p, which require formu- 
lating a renormalization-group transformation in terms 
of quenched probability distributions, we present the dis- 
tinct critical behaviors of the limiting cases of p = and 
p = l. 



A. Critical Properties at p = 

The d = 2, b = 2 Migdal-Kadanoff recursion rela- 
tions are exact |l7J on the p = lattice, and the 
renormalization-group transformation consists of deci- 
mating the two center sites in the cluster shown on the 
bottom right of Fig. ^ The renormalized Hamiltonian of 
the two remaining sites i' , j' is 



-PH' = J2 [J'si'Sj' +H' B {s v 

(i'j') 



i' 



J* = of the disordered phase and at low temperatures 
J > J c to the sink at J* = oo of the ordered phase. The 
critical behavior at T c is obtained from the eigenvalues 
of the recursion matrix at the critical fixed point, 



'2u CP 
2 + 2m u I , (17) 
1, 



where u = tanh2J c . This recursion matrix has eigenval- 
ues 2u = b VT , 2 + 2u = b VH , and 1, with eigenvalue expo- 
nents j/t = 0.747, yn = 1.879. Along the corresponding 
eigendirections are one thermal and two magnetic scaling 
fields: t = = hi = (2+ coth 2 J C )H B + H s , 

and hi 



1 


dJ' 


a.j' 


a.j' 


dJ 


dH B 


9H S 
dH B 




9H' B 


8H% 




dJ 


0H B 


dHs 
dH 1 s 


\ 


dH' s 


c)H' s 


OJ 


dH B 


dH s 



Jc ~ T c 

i?5, with linearized recursion relations t' 



b VT t, h[ = b yH hi, and h! 2 = h,2- Standard eigenvalue 
analysis at the fixed point yields the critical behaviors 
for the internal energy U = jj—^2(ij\ ( s i s j), the mag- 



netizations Mb, Ms, and the correlation length £: 



U-U c ~ \t\ 



i-, a = ^^ = -0.677, 
yr 



where the renormalized interaction constants are |4f 



M s ,M B ~\tf (t<0), /3=^-J^ = 0.162, 

yr 

£~\t\- v , v = — = 1.338. (18) 
yr 



J' = -hx{R ++ R—/Rl_) , 
H' B = ± In (R ++ /R— ) , H' S = H S , 
G' = 4G + - In (R ++ R—R 2 + _) , 



(14) 



with 



R++ = xy 2 z + x 1 z 1 , R = x 1 

R^ = yz + y~ 1 z~ 1 , x 



z + xy 2 z 1 



e 2J , y = e 2HB , z = e Hs 



(15) 



Here G is an additive constant per bond, equal to zero 
in the original Hamiltonian, but always generated by the 
transformation and necessary for the calculation of den- 
sities and response functions. From the transformation 
in Eqs. I|14f> . (|15fl we sec that, an initial Hamiltonian with 
only an Hs magnetic field term will invariably generate 
an Hb term upon renormalization. 

The subspace Hb — Hs — is up-down symmetric in 
spin space and closed under the transformation. Within 
this subspace, there is one unstable fixed point at 



Jr. = In 



- (l + (19 - 3\/33) 1/3 + (19 + 3V33) 1/3 



(16) 

corresponding to a temperature T c = 1/J C = 1.641. Un- 
der renormalization-group transformations, the system 
renormalizes at high temperatures J < J c to the sink at 



Mb and Ms have the same critical exponent /3, because 
the dominant magnetic scaling field h\ mixes Hb and 
Hs- Similarly, the susceptibility critical exponent is 7 = 
(2yH — d)/yx — 2.353. Approaching criticality in the 
ordered phase, all three susceptibilities one can define, 

Xbb - qh^, Xbs - \J — -qhj, and X ss - ohJ, have 
the critical behavior |i|~ 7 . The zero-field susceptibilities 
are infinite throughout the disordered phase. To recall 
this, we briefly review the calculation of thermodynamic 
densities and response functions by multiplications along 
the renormalization-group trajectory. 

Let K = (G, J, Hb, Hs) be the vector of interaction 
constants in the Hamiltonian, and K' = (G' , J' , H' B ,H' S ) 
the analoguous vector for the renormalized system. Cor- 
responding to each component K a of K is a thermo- 
dynamic density M a — ^- d ^ z , where Z is the par- 
tition function, and N a is a component of the vector 
N = (N nn , N nn , N nn , N n ). Thus, the density vector 
M = (1,U, Mb, Ms) is related to the density vector of 
the renormalized system M' by the conjugate recursion 
relations l47ll: 



M a = b 







M'T 0a 



Tf3a — 



N a dK a 



(19) 



An analogous recursion relation for response func- 
tions Xq/3 = \ / jf- §7^ has been derived by McKay and 
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Property 


p = 


< p < 0.494, 

(7 = 


0.494 < p < 1, 

(7 = 


<p < 1, 

< C7 < 1 


p=l, 

< a < 1 


<p < 1, 

(7 > 1 


T c 


1.641 


varies with p 
(see Fie [T2T 
reaches 3.592 
at p = 0.494 


varies with p 
fsee Fig IT2l- 
reaches 7.645 
at p = 1 


varies with a 
and p 
(see Fig. 


varies with a 

/qpp Pin- ISI- 

pec r ig. |0|l , 

reaches 3.485 
at cr = 1 


varies with a 
and p 
(see Fig. |SJ 


Vt 


0.747 
1.879 


vary with p 
(see Fig. [T4l 




1.585 


0.747 
1.879 


0.747 
1.879 


0.747 
1.879 


s 


£ -l/i/T 


£ -i/yr 


e^/Vm (t < 0) 
oo (t > 0) 


|C ^+/i( P ,^) 

1 1 


— -Ol(- In lil)"" 

1 1 


^ -i/ar 


C sing 






| i |-3/2g-2A/^ 










£ »T TJ2 ^' ' ' 




|t| «T 


M B , M s (t<0) 






e _A(2-3/ H )/vTt| 


^+^'" 




i "T 


XBB,XBS, , . 


\t\ y T 


^y/f — ^ 

t VT 


e A(2y H -2)/^\t\ 




N -^+ C 3(-l«| t |)-" 


\t\ »T 



TABLE I: Critical properties of the hierarchical-lattice network for all cases of long-range bond probabilities p and long-range 
bond decay exponents a. The functions fi(p,a,t) are corrections to the p = scaling behavior, which are significant at larger 
p. The factor A is defined below Eq. 1371 . and C\, C2, C3 are given as functions of yT,DH,cr in Eqs. (1621 . 



Bcrker kt 



v \ N X N^ , dK' x dK^ 
{;V N a N p Xx »dK a dK 



E 



(20) 



Using the density- response vector V = (I, U, Mb, Ms, 
Xbb, Xbs, Xss), Eqs. ljT§|) and (|2*U|l are combined into a 
single recursion relation, 



pa ■ 



(21) 



/3 



The extended recursion matrix W for the subspace Hb = 
H s = is 



(b d 






dG' 
8J 

dJ' 
dJ 










£>H' B 
9H B 








2dH^ 

8H' S 
Oil's 









e 2 C 
a 2 j' 






^dH B dH s 
u S 2 J' 
^dH B dH s 







2 e 2 g' 






( SH' B 
\8H B 



dH' B 8H' B 
^dH B dH s 
dH' B dH' s 
BH B dHs 



tr 1 

f* \dH s 
BH' B 9fl j 

\U " " " » (ffi 

(22) 

where ^ = y/N nil /N n . At a fixed point, V = V = V*, 
so that V* is the left eigenvector with eigenvalue b d of the 
extended recursion matrix evaluated at the fixed point, 
W*. To evaluate V for an initial system away from the 
fixed pont, Eq. I|21|l is iterated along the renormalization- 
group trajectory, 



(23) 



where V^™^ is evaluated in the system reached after the 
nth renormalization-group step, at which W' n ' is evalu- 
ated. When the total number of renormalization-group 
steps n is large enough so that the neighborhood of a 
fixed point is reached, ~V( n > ~ V* , so that V is evaluated 
to a desired accuracy, by adjusting n. 

From the recursion relations in Eqs. H14|) . l|15fl . the ex- 
tended recursion matrix W is 



w 



(A 2u 

2u 











2 + 2u ^ 

1 






\o 



iv 
-4m 2 





vEv 
VGu 2 







In 

3u 2 
2 






(2- 



■_ 

2 



J 



f2u) 2 V§u(l + u) 
2 + 2u 

1 

(24) 

1 + sech 2 J. At the sink of the 
disordered phase, u = 0, v = 2, and the left eigenvector 
of W* with eigenvalue b d is 



where u — tanh2J, v 



V* =(1, U = 0, Mb = 0, M s = 0, 

Xbb = 00, xbs = V6, Xss = 1) ■ (25) 

The matrix multiplication of Eq. (|23(l mixes xbb, Xbs, 
and Xss- Since xbb = 00 at the sink, all three sus- 
ceptibilities are infinite within the disordered phase. In 
contrast, at the sink of the ordered phase, u — 2, v = 1, 
and the two left eigenvectors of W* with eigenvalue b d 
are 

V* ± =(1, U=l, M B = ±2, M s = ±1, 

Xbb = 0, xbs = 0, xss = 0) . (26) 

Consequently, the susceptibilities from Eq. I|23|l are finite 
within the ordered phase, decreasing to zero as zero tem- 
perature is approached and increasing as \t\~ 7 as T c is 
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approached from below. The double value in Eq. (|26fl re- 
flects the first-order phase transition along the magnetic 
field direction. 

The infinite susceptibility in the disordered phase is di- 
rectly related to the presence of sites with arbitrarily high 
degree numbers in the scale-free network, because these 
sites feel a very large applied field, channeled through 
their many neighbors. Except for this feature, the critical 
behavior for the p = case is similar to that of a regular 
lattice, which is unsurprising since the Migdal-Kadanoff 
recursion relations that are exact on the hierarchical lat- 
tice can be derived from a bond-moving approximation 
applied to the square lattice. 

The p — results are in Fig. 03 where the specific 
heat, magnetizations, and zero-field susceptibilities are 
plotted as a function of temperature. Since the specific 
heat exponent is a — —0.677, the specific heat has a finite 
cusp singularity at T c . 



B. Critical Properties at p = 1 

For the p = 1 lattice, the renormalization-group trans- 
formation consists of decimating the two center sites in 
each connected cluster of the type shown on the top 
right of Fig. ^ The Hamiltonian now includes long-range 
bonds, Eq. (|11J) . and the transformation is a mapping of 
the Hamiltonian —(3H.(J, Hb, Hs, {K q }, G) onto a renor- 
malized Hamiltonian -j3'H'(J', H' B ,H' S , {K' q }, G 1 )). The 
recursion relations are 



J' = \\n(R, 
H> B = \\n{R+ 
G' = AG + \ In (R ++ R—R 2 + _) , 



-R — /R\. 

-/R-) , 



9 = 1,2,. 



(27) 



where R ++ , R , and i? H are as given in Eq. (|T5)) . 

Long-range bonds as well as nearest-neighbor bonds 
now contribute to the internal energy U, 



U = 



N nn U nn 



(28) 



where 



C/nn = 



1 d 



hxZ. 



U\d,q — 



ld,q 



E 



N nn dJ 
_ _l d_ 

' N M ,g dK C 



InZ. 



(29) 



Here N\d,q = Ar q N nn is the number of long-range bonds 

with rriij = q. Since K' q does not depend on J, Hb-, or 

Hs, the thermodynamic densities and response functions 

in V = {l,U nn ,M B ,M s ,XBB,XBS,Xss) still obey the 

< — > 

recursion relation in Eq. I|21|) with a matrix W of the 
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FIG. 5: Specific heat, magnetizations, and zero-field magnetic 
susceptibilities for p = 0, as functions of temperature 1/J. 
The dotted vertical line marks the critical temperature T c = 
1.641. As insets to the magnetizations and susceptibilities, we 
show InM and lnx with respect to ln|t|, where t = ^ " < 
0. The linear behavior in the insets agrees with the power- 
law predictions of M B , M s ~ |t| 0162 and xbb, Xbs, Xss ~ 

| t |-a.36S > 



same form as in Eq. (|22|l . The densities l/id ig , on the 
other hand, have the recursion relation 



Uua = b- d U' D 



N nD dJ' = 
'JVidiMi 



b~ d Ui 



Id, g-1 " 



Id, q 



dK a 



= UL, q -i (9>1)- 
(30) 



Thus U\d, q — Unn , where Unn is the nearest-neighor den- 
sity U nn in the system reached after q renormalization- 
group transformations. Thus all the long-range bond 
densities Uid, q are found by evaluating U nn along the 
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FIG. 6: Three possible behaviors of the renormalization-group flows of the p = 1 network with uniform long-range bonds. The 
curve in each diagram is the recursion J' (J) from Eq. i>! . with the straight line J' — J also drawn for reference. Intersections 
of the curve with the straight line are fixed points. The flows are given by the staggered line, with successive values of J' 
corresponding to where the staggered line touches the curve. Only the dotted fixed points are physically accessible. The inset 
on the right shows the continuous line of fixed points J* (Jo) as a function of Jq. 



renormalization-group trajectory. Eq. (|28|) can be rewrit- 
ten as 



2.9 

9=1 



(31) 



where we have also used N nn +J2^Li Md, 9 = f Ann- From 
Eq. lpTT|l and the recursion relation for £/„„, the leading 
singularity in U nn is also the leading singularity in U . It 
is sufficient to calculate the singular behavior of U n n to 
obtain the critical properties of U and of the specific heat 
C. 



1. Long-distance bonds with uniform interaction strengths 

We first consider the case with no distance depen- 
dence in the strengths of the long-range bonds, a — 
0. Here K q = Jq for all q and after any number of 
renormalization-group transformations, where Jo is the 
value of J in the original system. The recursion relation 
for J in the closed subspace Hb — H$ — is 



J' = J + ln(cosh2J) . 



(32) 



There are three types of behavior possible for the 
renormalization-group flows, as illustrated in Fig. For 
Jo greater than a critical value J c (Fig. EJa)), the flows 
go to the ordered phase sink J* — oo. For J < J c 
(Fig. EIb,c)) the flows go to a continuous line of fixed 
points J* (Jo), with a distinct fixed point for each start- 
ing interaction Jo. When Jo = J c exactly, the J' (J) curve 
touches tangentially the straight line J' = J at J*(J C ), 
as shown in Fig. E{b). This fact allows us to solve for 
J*(J C ) and J c exactly: 



J*(J C ) = -ln3, J c = ln 



3 3/4 



(33) 



Thus the system is conventionally ordered below the 
critical temperature T c = 1/J C = 7.645. To under- 
stand the novel high-temperature phase above T c , we 
look at the recursion matrix W* evaluated along the 
line of fixed points, J* (Jo) for Jo < J c . The form of 
the matrix is as in Eq. J2U, with u = tanh2J*(Jo) and 
v = 1 + sech 2 2 J* ( Jo). Since J* (Jo) has the maximum 
value of (ln3)/4 = 0.275 for Jo = J c and tends to zero 
as J increases, < u < 1/2, 7/4 < v < 2. The left 

eigenvector of W* with eigenvalue b d is 



=(1, u n 



2-u 



M B = 0, M s = 0, 



Xbb = oo, xbs = oo, XSS = oo) ■ (34) 

It follows that, in the high-temperature phase, Mb = 
Ms = and that the susceptibilities xbb, Xbs, XSS are 
infinite. Because the renormalization-group flows go to 
a line of fixed points ending at the critical point J*(J C ), 
the correlation length is infinite throughout the phase 
and the correlations have power-law decay, characteris- 
tics which are typically seen just at T = T c . (In contrast, 
the low-temperature ordered phase has the usual expo- 
nential decay of correlations.) This type of behavior, with 
a transition between phases with finite and infinite corre- 
lation lengths, was first seen in the Berezinskii-Kosterlitz- 
Thouless phase transition 29, 30], though with an impor- 
tant difference: There the algebraic order was in the low- 
temperature phase, while here it is the high-temperature 
phase that has this feature. 

We now turn to the critical behavior of the system in 
the ordered phase, as T — » T c from below. For small 
negative t = (T — T c )/T c = (J c — Jq)/ Jq, we have Jo 
,/, + S, where S = J c \t\- As can be seen from Fig. (6a), a 
renormalization-group flow starting at Jo spends a large 
number of iterations in the vicinity of J*(J C ) = (ln3)/4, 
before escaping to the ordered phase sink at J* = oo. If 
no is the number of iterations initially required to get J 
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close to J*(J C ) and n* is the number of iterations where 
J rs J*(J C ), then as S — > 0, n remains constant, while 
n* — > oo. The dependence of n* on (5 (and hence on 
determines the critical singularities. For a typical 
critical point, n* ~ (ln<5)/(?/T In 6). However, in our case, 
at J*(J C ) the eigenvalue exponent j/t = 0, and it turns 
out that n* ~ S^ 1 / 2 . We show this as follows: After no 
iterations, the flow is at J near J*(J C ), with J < J*(J C ). 
It then takes n*/2 iterations to get J almost exactly at 
J*(J C ), and another n*/2 iterations to get J a significant 
distance away from J*(J C ), namely to J - J c ~ 0(1). 
Considering the latter half of this flow, we expand the 
recursion relation for J, Eq. (|32|l . around J*(J C ), 



J' - J* (Jo) = 5+ (J -J* ( J c )) + 2(J - J* (J c ))2 
- (J — J*(J C )) 3 + • • • . 



(35) 



Starting with J = J*(J C ), from Eq. I|35j) . we ob- 
tain series expressions for jW, the interaction after i 
renormalization-group steps: 

J (1) - J*(J C ) = 5 

J*(J c ) = 2S+^5 2 - 5 3 + --- 



j(n) _ j* j =n g + ~(n- l)n(2n - 1)<5 2 

+ — (n - 3)(n - l)n(24n 2 - 29n + 2)<S 3 + • • • 
80 

(36) 

For n -C <5 _1 / 2 , the first term in the series for J^™- 1 — 
J*(J C ) is dominant, and the distance increases very 
slowly as j( n ' — J* ( J c ) ~ n5. For large n, the kth term in 
the series ~ n 2k ~ l 5 k . Thus, when n is of the order <5~ 1//2 , 
_ J*[J C ) begins to increase significantly. From this 
we can deduce that n* scales like 5~ 1 ^ 2 . 

We can now proceed to find the critical behaviors for 
the correlation length, thermodynamic densities, and re- 
sponse functions. By iterating the recursion relation for 
the correlation length, £ = £ = b n * +n ° £ (n+n °\ where 
is the correlation length after n renormalization- 
group steps. The singularity in £ as 5 — > comes from 
the 6™ factor, 



£~b n 



C In 2 

e 



A 



(37) 



where n* sa C<5 for some constant C, and A = 

c/VT c . 

From Eqs. (|22l) . 1231) . we extract the critical behaviors 
of the internal energy, magnetizations, and susceptibil- 
ities: The nearest-neighbor contribution to the internal 
energy U nn transforms as 



Uu 



k-d 



dG 
~~dJ 



k-d 



B V 
U' — 

nn dj 



(38) 



Since dG' /dJ is analytic, the singularity of U nn must 
reside in tAm sing = b~ d U' nYy dJ'/dJ. Iterating over n + n* 



renormalization-group steps, 



n +n* 



l— n*d 



i=l 
n 



dJ' 
~dJ 



h -n dTj(n a +n') TT UJ 

11 Q J 



J=JW 



J=JW 



(39) 



where we have used the fact that OJ'/dJ ~ 1 for the 
n* iterations during which jW ~ J*(J C ). After n + n* 
iterations the system has flowed away from criticality. 
The singular dependence comes from the b~ n d factor, 



U n 



-n d 



Thus the singular part of the specific heat is 



The magnetizations Mb and Ms recur as 



(M B) M s )=6- d (M^,M^) 







1 



(40) 



(41) 



(42) 



where u — tanh 2 J. Iterating over no + n* 
renormalization-group steps, 

{M B , M s ) 

~ 6-("«+"*) d ((M^ 0+n * ) ,M|" 0+ "* ) ) • v) 6™*^v R. 

(43) 

Here 6 yff = 2 + 2 tanh 2 J*(J C ) = 3 is the largest eigen- 
value of the 2x2 derivative matrix in Eq. I|42|l evaluated 
at J*(J C ), v is the corresponding normalized (to unity) 
eigenvector, and R is the product of the derivative ma- 
trices of the first no iterations. The singular behavior 
comes from the factor b~ n ( d ~y H \ 

M B , M s ~ b- n ' (d - yH) ~ e _< . (44) 

Since yn — 1.585, the magnetizations decrease exponen- 
tially to zero as |t| — * 0. 

The susceptibilities xbb, Xbs, an d Xss recur as 



(xbb,Xbs,Xss) 



— b {x'bbiX'bS'X'ss) 



'(2 + 2k) 2 V6u(1 + u) 2|i 



2 + 2u 



+b d (G\K n ) ( jj 2 . 



(45) 



where v = 1 + sech 2 J. Since there is no singular behav- 
ior in G" and U' — s- as Itl — ^ 0. only the first term 

nn sing 11 

in Eq. I|45fl contributes to the divergent singularity of the 
susceptibilities. Iterating over no + n* steps, 

(XBB,XBS,X55)si„g~&~ (,l0+n * )d 



( ™ 0+n *\xBl + "^Xss + "* ) ) • v ) 6 2 "* yif v R, (46) 
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where b 2yH = (2 + 2tanh2J*(J c )) 2 is the largest eigen- 
value of the 3x3 derivative matrix in Eq. (|45|) evaluated 
at J*(J C ), v the corresponding normalized eigenvector, 
and R the product of the derivative matrices for the first 
no steps. The singularity in the susceptibilities is given 
by 



0.3 



XBB 



sing " 



XBS S 



XSS S 



(2y H -d)A 



(47) 

We illustrate these results in Fig. plotting the spe- 
cific heat, magnetizations, and zero-field susceptibilities 
as a function of temperature. The essential singularity 
in the specific heat (Eq. is invisible in the plot, the 

function and all its derivatives being continuous at T c , 
with the rounded analytic peak occurring in the phase 
opposite to the algebraic phase, namely in the ordered 
phase at lower temperature. This behavior of the spe- 
cific heat also occurs in the XY model undergoing a 
Berezinskii-Kosterlitz-Thouless phase transition, as seen 
in Fig. 5 of Ref.[48|. In the latter case, opposite to the 
algebraic phase, the phase in which the rounded analytic 
peak occurs is the disordered phase at higher tempera- 
ture. In the XY model, the physical meaning of the high- 
temperature rounded peak is the onset of short-range or- 
der within the disordered phase. In our current system, 
the physical meaning of the low-temperature rounded 
peak is the saturation of long-range order that occurs 
unusually away from criticality, due to the essential crit- 
ical singularity of the magnetization, which corresponds 
to a critical exponent [3 — oo and the unusual flat onset 
of the magnetization, as seen in Fig. [7| 



2. Long-distance bonds with decaying interaction strengths 




4 6 8 

Temperature 1/J 
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For a > 0, the long-range bond strengths K q = J q a 
and thus at the nth renormalization-group step — 
Jo(n+ l) _cr . The interaction strength j( n ) in the closed 
subspace Hb — H$ — is given by the recursion relation 

J {n) = J n- a + ln(cosh 2J {n ~V) , (48) 

where J^ -* = Jq. The critical temperature T c now de- 
pends on the exponent a, as shown in top curve of Fig- EI 
having the maximum value of T c — 7.645 at a = and 
decreasing with increasing a (to T c — 2.744 at a = oo, 
where the system reduces to a nearest-neighbor, next- 
nearest- neighbor model) . 

When the number of renormalization-group steps n — > 
oo, the Jon~ a term in Eq. (|48|l goes to zero, so that the 
fixed points of the renormalization-group transformation 
are those of the p = case analyzed in Sec. III. A. Thus 
for temperatures close enough to T c , satisfying \t\ <C r for 
some crossover value r, we expect to observe the p = 
critical behavior. However, the width r of the critical 
region varies with a, becoming extremely narrow as a — > 
0. For a thermodynamic quantity scaling as \t\ x inside 



FIG. 7: Specific heat, magnetizations, and magnetic suscep- 
tibilities for p — 1 long-range bonds with uniform interaction 
strengths (a = 0) , as a function of temperature 1/J. The dot- 
ted vertical line marks the critical temperature T c — 7.645. 
The insets to the magnetization and susceptibility graphs 
show In | lnM| and ln(lnx) versus In \ t\, where t = ~ Te < 0. 
The linear behavior in the insets agrees with the exponential 

scaling predictions of Mb, Ms ~ e~ c ^^^ and Xbb, Xbs, 

Xss ~ e D// V^ with positive constants C and D . 



the critical region [x being one of the p = exponents), 
the general scaling behavior for small \t\ not necessarily 
in the critical region is !*!*+/»(*), where \f x {t)\ < |x| 
when |t| -C t, and the form of f x {t) may depend on a. 
In the following, we derive the leading order contribution 
to f x (t) for the various physical properties of the system, 
also determining the size of the critical region r. 

If the system is at its critical temperature, Jo = J c , 
the interaction strength under repeated renormalization- 
group iterations, j' n J for n — > oo, goes to the p = crit- 
ical fixed point, which we will label J c o and whose value 
is given by Eq. I|16|) . Let us denote this renormalization- 
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FIG. 8: Critical temperature T c = 1/J C for various p as a 
function of the exponent a describing the decay of the long- 
range bond interaction strengths. The curves for < p < 1 
were calculated using the techniques in Section III.C. 



group flow as Jc™\ so that Jc = J c and lim^oo — 
J c o- Now if we start instead at a temperature very close 
to critical, J = J c —J c t for small \t\, j(") stays near j£ n ' 
for a large number of iterations n* , before veering off to 
either the ordered or disordered sink. The dependence 
of n* on \t\ is the key to the crossover behavior of the 
system. The difference — ji n ^ satisfies the recursion 
relation 



j(n+l) _ j(n+l) = 61W (n)(j(n) _ j(n)j ) (49) 



where 



dj(n+l) 



= 2 tanh 2 jj n) 



Iterating Eq. iTTOl) . 

j(n+l) _ J( (n+1) = b E n k=0 VT(k)(j o _ 
= &E^ =0 i/TW(_j c t). 

Since J(»*) - J c (n *> - 0(1), 

In If I 



(50) 



jSJ=oiw(*) ^ an d 



k=0 



In 6 



(51) 



(52) 



In order to find n* , we need to determine j/y (n) . From the 



fact that lim 



n — >oo ° c 



(n) 



J c o and the recursion relation 
in Eq. (|48|l . we consider for Jc the large n form of 

4") = J c0 - + • • • . (53) 

Substitution into Eq. lj4"5|) yields 

Jo 



B = 



2 tanh 2 J c , 



1 



(54) 



Eqs. (|53|l . (|54|) can also be obtained by expanding the 
recursion relation around J c o, 

4 n) - J c0 = J n- a + 2 tanh 2 J c0 ( jj"" 1 ) - J c0 ) , (55) 



and summing the series derived from iterating Eq. i|55[l . 
Substituting into Eq. lj5U|) . 



2/t(«) =yro-Cn a + 



(56) 



where ?/to = ln(2tanh2J c o)/ln& = 0.747 is the p = 
thermal eigenvalue exponent and C = 1.498Jo- For use 
below, we also deduce the magnetic exponents y#(n), 



l)VH(n) 



dH 



(n+l) 
B 



H 



(n) 



(») rrCn). 



(57) 



= 2 + 2tanh2J,(™), 

(n) = J/ffo - BrT a + ■ ■ 

where y H o = ln(2 + 2tanh2J c0 )/lnfe = 1.879 is the p = 
magnetic eigenvalue exponent and D = 0.683 Jo- 

From Eq. j55} , we evaluate X)fc=oy T (^) f° r l ar g e n *j 




CC(<T) : 

For a > 1, the n*yTo term is clearly dominant for large 
n*, so that, from Eq. 

1 In 1*1 



< a < 1 
a = 1, 
CT > 1. 



(*>1) 



(59) 



" — ; — r = "0 

Uto mo 

This expression for n* leads to the same critical expo- 
nents we found in the p = case. On the other hand, for 
the slow decay of < a < 1, Eq. I|52|l becomes 

Cn' 1 -" In |t| 
n yro : - — r~ T • ( 60 ) 



1 



In b 



Writing n* = Uq + 5n, the leading order contribution to 
Sn is found. 



Cn* 



(1 - a)y T o 



(0 < a < 1) 



(61) 



This expression for n* when < a < 1 yields the 
leading-order corrections to p = in the critical behav- 
iors of the correlation length, internal energy, specific 
heat, magnetizations, and susceptibilities: 



Using 
Csing 

M B , M s 

XBB, XBSi XSS 



b n ~ |*| VT0 

5-™* d +Efc =0 s/T(fe) 



f*=-(- 



l°|t 
In b 



"TO (l-a)y 



ac ( In |t \- 



d — 2y TQ 



TO 
dC 



1*1 
1*1 



• TO 



TO 



1*1 

l ) -n*d+2'£k=a VH(k) 



-2BJJQ (d-2y gn )C + 2 HT0 -D 



»T0 



0--")v% 



(62) 
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All the critical behavior expressions in Eqs. (|62|l have 



the form |if +<1 



where x is the appro- 
priate p = exponent and E is a non- universal (i.e., Jo- 
dependent) constant ~ 0(1). For temperatures \t\ < r, 
the leading-order correction term in the exponent should 
be negligible, 



1 



(i - °)vt7 



ln|f| 
In 6 



< 



(63) 



for some small quantity e, giving an estimate for r as 

(7 — > 0, 



ft -(e(l-cr)l&-")- 



(64) 



With decreasing a the critical region r becomes rapidly 
infinitesimal. For example, with a — 1/2 and e = 10 , 

T » 10~ 289 . 

The above corrections to critical behavior are illus- 
trated in Fig. El where we plot numerically calculated 
effective exponents lnMss/ln|t| and In xbb I In \t\ as a 
function of \t\ for several values of a. It is clear that 
for a > 1, the effective exponents quickly converge to 
the horizontal lines showing the actual asymptotic expo- 
nents. The convergence when a < 1 is much slower, due 

E ? I \n\t\ y" 

to the \t\ ''""'•to" " correction to asymptotic uni- 
versal critical behavior. In Fig. llUl we explicitly show for 
the case a = 0.6 the magnetizations and susceptibilities 
asymptotically approaching the scaling forms of Eq. I|62|l 
for small \t\. 



C. Critical Properties of the System with 
Long-Range Quenched Randomness, < p < 1 

1. Exact renormalization- group transformation for 
quenched probability distributions 

When < p < 1 , there is long-range quenched random- 
ness in the network, and the renormalized system will 
have an inhomogenous distribution of all interaction con- 
stants. The renormalization-group transformation needs 
be formulated in terms of quenched probability distribu- 
tions |4Jj . First consider the decimation transformation 
effected on the cluster of Fig. ^2 with nonuniform in- 
teraction constants. The recursion relations for J'i'j>, 
H' Bi ,-,, H' s , and G' v y are the locally differentiated ver- 
sions of Eq. J23, 

J'ii'j') = J'd'hf) + J'(i'k 2 f) + Kiti'f) , 

j'(i>k l j') = \ln{R ++ R--/Rl_) {i , kxjl) , 

H' B (i'f) = H'xd'kd') + H' B (j%f), H' s (i') = H s (i') , 

H B (i'k 1 j') = lln(R ++ /R—) (i , kijl) , 




10 -70 10 -60 10 -50 



1*1 



FIG. 9: The calculated effective exponent of the magnetiza- 
tion Mb and magnetic susceptibility Xbb, as a function of \t\ 
for t — T ~ Tc < and p = 1. Curves for several values of a, 
the exponent for the decay of the long-range bond strengths, 
are shown. The horizontal dashed line in the upper graph 
corresponds to the actual critical exponent for the suscepti- 
bility, —7 = — (2yH0 — d)/j/T0 = —2.353, while the dashed line 
in the lower graph corresponds to the actual magnetization 
exponent, f3 — (d — yro)/yTO = 0.162. 



G'^k 1] ') = - A \n{R ++ R^R\_) [i , kijl) , 

K' q = K g+1 , g =l,2,..., (65) 

where R++, R — , and R-\ — along path (i'kj') are given 
by the locally differentiated versions of Eq. lfT5l> , 

R++ = x x x 2 y\y 2 z + x^x^z^ 1 , 



R = x 1 x 2 z + £1X22/1 V 2 z 

R+- = xix^yiz + x 2 y 2 z -1 , 
Xl =e J ^' k \x 2 = e J ^'K 

m = e 2H B (i'fc) ; y2 = e 2ff B(fc /) j z = e H s (k)_ (6g) 

If there is no long-range bond connecting i' and j' , the 
equations above hold with K\(i! j') = 0. We shall work 
in the closed subspace H B (ij) — Hs{i) = for all 
where the recursion relation for J'(i'j') is a function 



J'(i'j')=R({J(ij)};K 1 (i>j')) 



(67) 



G , (i , j , ) = G , (i , k 1 j , ) + G , (i , k 2 j'), 



with {J(y)} = {J(i'fci),J(Ai/). J ( i/fe 2), J(k 2 f)} being 
the set of interaction constants in the cluster, and R given 
in Eqs. (jfJ5|) and lj6U|) . 
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FIG. 10: The effective exponents of the magnetizations Mb 
and Ms, and susceptibilities Xbb, Xbs, and xss, for p = 1, 
a — 0.6, plotted as a function of |i| for t — T ~ Tc < 0. 
The dotted curves in the figures are the first-correction-to- 
scaling predictions for the magnetization and susceptibility 
from Eq. 1621 . The horizontal dashed line in the upper graph 
corresponds to the p — critical exponent for the susceptibil- 
ity, —7 = — (2?/£ro — d)/yTo = —2.353, while the dashed line 
in the lower graph corresponds to the p = magnetization 
exponent, f3 = (d- yro)/yTo = 0.162. 



If the interaction constants J(ij) have a quenched 
probability distribution V(J(ij)), and the long-range 
bonds K q (ij) have a quenched probability distribution 
Q ig) (K q (ij)), the distribution V {n) {J' v f ) for the rescaled 
system after n renormalization-group transformations is 
given by the convolution 

v^(j'a'j')) = 

[II dJ(ij)v^ ( j(y))l m c (B - 15 (#1 (i'i')) 



5 (J'i 



(68) 



where the product runs over the nearest-neighbor bonds 
ij in the cluster between i' and j' . The long-range bond 
distribution Q( n \K[ (i'j 1 )) after n renormalization-group 
transformations is 

Q^(K[(i'j')) = P 6 {K'S'j') - J (n+ 1)- CT ) 

+ (l-p)<f(^ (t'j')). (69) 

The convolution in Eq. (|68|l is implemented numer- 
ically, with the probability distribution V^ n \Jij) rep- 




FIG. 11: Cluster with quenched randomness on which the 
decimation transformation of Eq. 16511 is applied. 



resented by histograms, each histogram consisting of a 
bond strength and its associated probability. The initial 
distribution V (0 \Jij) is a single histogram at Jo with 
probability 1. Since Eq. (|68|l is a convolution of five prob- 
ability distributions, computational storage limits can be 
used most effectively by factoring it into an equivalent 
series of three pair wise convolutions, each of which in- 
volves only two distributions convoluted with an appro- 
priate R function. Two types of pairwise convolutions 
are required, a "bond-moving" convolution with 



Rbm{J{hh),J{i232)) = J(kji) + Jihh) ■ 



and a decimation convolution with 



(70) 



In 



cosh(J(iiji) + J {1232)) 



2 Vcosh(J(iiji) - J(hh)), 

(71) 

Starting with the probability distribution V^ n 1 \ the 
following series of pairwise convolutions gives the total 
convolution of Eq. (|68|) : (i) a decimation convolution of 
pi"- 1 ) with itself, yielding Va\ (ii) a bond- moving convo- 
lution of V a with itself, yielding Vb'i (iii) a bond- moving 
convolution of Vb with Q'™ -1 -', yielding the final result 

Because the number of histograms representing the 
probability distribution increases rapidly with each 
renormalization-group step, we use a binning proce- 
dure |50j: before every pairwise convolution, the his- 
tograms are placed on a grid, and all histograms falling 
into the same grid cell are combined into a single his- 
togram in such a way that the average and the standard 
deviation of the probability distribution are preserved. 
Histograms falling outside the grid, representing a neg- 
ligible part of the total probability, are similarly com- 
bined into a single histogram. Any histogram within a 
small neighborhood of a cell boundary is proportionately 
shared between the adjacent cells. After the convolu- 
tion, the original number of histograms is reattained. For 
the results presented below we used 562,500 bins, requir- 
ing the calculation of 562,500 local renormalization-group 
transformations at every iteration. 

For the thermodynamic densites M a , given by 



d\nZ 
\ n ^ dKjrj) 



(72) 



the chain rule yields conjugate recursion relations for the 
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quenched random system, 



cr = 



h -d V — V dlnZ V 



, N 'p fp 9K' (i'f) ^ N a dK a (ij) ' 



where the rightmost sum runs over nearest-neighbor 
bonds ij in the cluster between sites i' and j'. As an 
approximation, this sum is replaced by its average value, 
so that 



Here the overbar denotes averaging over the probabil- 
ity distributions of the interaction constants in the clus- 
ter shown in Fig. 1111 Using the recursion relations in 
Eq. iJBSJ, in the subspace H Bi j = H s = 0, 




0.2 0.4 0.6 0.8 

Probability of long-distance bonds p 

FIG. 12: Critical temperature T c as a function of the proba- 
bility of long-range bonds p, plotted for several values of the 
exponent a characterizing the decay of the long-range bond 
strengths. 



(a X^i'j' dG'ji'j') 
I * L^ij dJ(ij) 

n V i/j ' dJ '( l 'j") 
u l^ij dJ(ij) 







\o 

/4 2u 
2u 







S-^i'j' dH' B (i'j') s-ii'j' dH' B (i'j') 

Z-^ij dH B (i 3 ) dHs(i) 

„ sr^i' 9H' s (i') 

u l^i dHs(i) 




(75) 



\0 



where 



u = -(tanh(J(i'fci) + J(hj')) 

+ tanh(J(« / fe 2 ) + J(k 2 j'))), 
sinh(J(i'fci) + J(fcu') + J(i'k 2 ) + J(k 2 j')) 



w = 



2cosh(J(i'fci) + J(fcij')) cosh( J{i'k 2 ) + J(k 2 j')) ' 

(76) 

For a fixed probability distribution of the 
renormalization-group transformation (e.g., Fig. 113(1 . the 
thermal and magnetic eigenvalues exponents tjt and yn 
are obtained as 



ij 
i'f 



dJ'ji'j') 
9J(ij) 



dH B (ij) 



2u, 



= 2 + 271. 



(77) 



2. Results 

The quenched random system critical temperatures 
T c (p) — J^ 1 {p) are shown as a function of p, in Fig. IT21 
for several values of the decay exponent a. For any a > 0, 
in renormalization-group trajectories starting near T c , 
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FIG. 13: Histograms of the unstable critical fixed probability 
distributions for p = 0.2, a = (left column) and p = 0.7, 
a — (right column). The bottom panels show the actual 
histograms (numbering 1,128,002 in each case), while in the 
top panels the histograms are combined in order to clearly see 
the outlines of the probability distributions. 



the probability distribution V(J(ij)) spends many itera- 
tions in the vicinity of the unstable critical fixed distri- 
bution which is a delta function at J c {p = 0), the p = 
critical interaction strength given by Eq. 1(16(1 . Similarly 
to the results of the p = 1 case given above, when a > 0, 
the critical behavior for all p is that of p = 0, though 
with a rapidly decreasing critical region as p — > 1 and 
a 0. 

On the other hand, for a — 0, a variety of critical be- 
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FIG. 14: The thermal and magnetic eigenvalues yr and yn, 
and the corresponding specific heat and magnetic exponents 
a and /3, as a function of p, for a = 0. The probability p = 
0.494, marking the onset of exponentiated power-law critical 
behavior, is shown with a dotted line. 



haviors occurs asp ranges from to 1. The unstable crit- 
ical fixed distribution has a non-trivial structure which 
depends on p, two examples of which are shown in Fig. 1131 
The eigenvalue exponents j/t and yn from the critical 
fixed distributions change continuously with p (Fig, 114}. 
and with them the critical exponents characterizing the 
phase transition. As p is increased from zero, both i/t 
and yn decrease from their p — values, attaining their 
p = 1 values of yr — and yn = In 3/ In 2 at p = 0.494. 
Thus, the system has two distinct regimes of criticali- 
ties. For p < 0.494 the critical behavior is described by 
power laws with exponents v = 1/yr, & = (2yr — d)/yr, 
f3 = (d~y H )/y T , and 7 = (2y H ~d)/y T - As y T -> with 
p — > 0.494, the exponents blow up as v — > 00, a — > —00, 
(3 — ► 00, and 7 — * 00. For p > 0.494 the critical behavior 
is that of the p = 1, a = case given above, with expo- 
nentiated power laws of the thermodynamic quantities, 
and the high-temperature phase has infinite correlation 
length. The onset of exponentiated power-law critical be- 
havior at p = 0.494, due to the influence of the long-range 
bonds, in fact corresponds to a change in the geometri- 
cal features of the network. As we have noted in Fig. 0] 
f° r P ^ 0.5 the average path length £ has a small world 
character, £ ~ In N n , while for smaller p it increases more 
rapidly like Nn , as in a regular lattice. 

The spectrum of critical behaviors for varying p at 
a = is illustrated in Figs. El an d El for the specific 
heat and magnetization versus temperature for different 
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FIG. 15: Specific heat calculated for various probabilities 
p, plotted with respect to the normalized temperature (T — 
T c )/T c , with a — 0. The vertical dotted line corresponds to 
the critical temperature T — T c . The inset shows a close-up 
of the specific heat near T c for p — 0.06, showing both the 
infinite-slope singularity at T = T c , and the analytic peak 
that appears for T < T c when p > 0.053. 




FIG. 16: Magnetization calculated for various probabilities p, 
plotted with respect to the normalized temperature variable 
(T c - T)/T c , with o = 0. For p > 0.363, note the unusual flat 
onset at T c . 



values of p. With increasing p from p — 0, the low- 
temperature analytic peak, due to the saturation of long- 
range order, as mentioned above, appears atp« 0.053, 
as the low-temperature amplitude of the critical cusp 
changes sign, and shifts to lower temperatures as p fur- 
ther increases. At the critical-point singularity, with in- 
creasing p from p — 0, the specific heat exponent a con- 
tinuously decreases from its p = value of —0.677: The 
cusp disappears at p = 0.105 as a crosses —1, so that 
the specific heat acquires a continuous slope at critical- 
ity, but all higher derivatives remain divergent. The sec- 
ond derivative at criticality also becomes continuous, all 
higher derivatives remaining divergent, at p = 0.249 as 
a crosses —2. Thus, as a crosses the consecutive nega- 
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tive integers at at p = 0.105, 0.249, 0.312, 0.349, . . ., the 
divergence begins at a higher derivative, until the accu- 
mulation point at p — 0.494, where a reaches — oo, and 
the essential singularity occurs for the higher values of p. 

In the magnetization, with increasing p fromp = 0, the 
critical exponent (3 continuously increases from its p — 
value of 0.162. Thus, the slope at criticality changes 
from infinity to zero at p = 0.363 as (3 crosses 1, but all 
higher derivatives of the magnetization remain divergent. 
The second derivative at criticality also becomes zero, all 
higher derivatives remaining divergent, at p — 0.424 as (3 
crosses 2. At each crossing of a positive integer by (3, at 
p = 0.363, 0.424, 0.446, 0.457, . . ., the zeros extend to one 
higher derivative and the divergence begins at one higher 
derivative, until the accumulation point at p = 0.494, 
where (3 reaches oo, and the essential singularity occurs 
for the higher values of p. 



IV. CONCLUSIONS 

In summary, we have introduced a scale-free 
hierarchical-lattice network model exhibiting a range of 
geometric and thermodynamic properties as we vary the 
probability p of the long-distance bonds. When p = 0, 
our network is unclustered and the average shortest- 
path length t scales like a power-law in the number of 
sites, I ~ N 1 / 2 , resembling in these respects a regu- 
lar lattice. This resemblance also holds for the crit- 
ical behavior of the ferromagnetic Ising model on the 
p = network, which shows typical power-law singu- 
larities in specific heat, magnetization, and susceptibil- 
ity. For small concentrations of long-distance bonds, 
p < 0.494, this picture does not change radically: the 
clustering coefficient increases nearly linearly with p, the 
average shortest-path length continues to have power-law 
scaling with lattice size, and the critical exponents of the 
Ising model vary continuously with p. As we approach 
p = 0.494, however, the magnitudes of these exponents 
blow up, and we have an unexpected crossover to a com- 
pletely different regime of critical behavior for p > 0.494. 
We find a highly unusual infinite-order phase transi- 
tion, an inverted Berezinskii-Kosterlitz-Thouless singu- 
larity between a low-temperature phase with nonzero 
magnetization and finite correlation length, and a high- 
temperature phase with zero magnetization but infinite 
correlation length and power-law decay of correlations 
throughout the phase. This slow decay of correlations 
in the disordered phase is a direct consequence of the 
underlying lattice topology, since large enough concen- 
trations of long-distance bonds significantly reduce the 
shortest-path length between any pair of sites. Indeed 
for p > 0.494 the network shows the small-world effect, 
with the average shortest-path scaling logarithmically as 
€~lniV. 

In determining these geometric features and critical 
behaviors of our model, we were aided by the deter- 
ministic, hierarchical nature of the network construction. 



This allowed us to derive analytic expressions for many 
of the network characteristics: the degree distribution 
and the clustering coefficient at all p, and the average 
shortest-path length for p = and p = 1. In addi- 
tion, we were able to formulate exact renormalization- 
group transformations for the Ising model on the net- 
work, even with a quenched random distribution of the 
long-distance bonds. The present model was designed 
to incorporate just some of the distinctive properties of 
real-world networks: scale-free degree distribution, high 
clustering coefficient, and small- world effect. But such 
hierarchical-lattice models could also be tailored to cap- 
ture other empirical properties, like the modular, com- 
munity structure of networks. Using techniques similar 
to the ones we applied to the Ising model, one could de- 
velop exact renormalization-group approaches to other 
interesting statistical physics systems, for example per- 
colation models relevant to epidemic spreading. Finally, 
non-equilibrium study [5lJ of our model should yield in- 
teresting new results. 
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APPENDIX A: DERIVATION OF NETWORK 
CHARACTERISTICS 

1. Average clustering coefficient C m 

Consider a site in the infinite lattice with k nn = 2 m , 
feld = 2 m — 2, and m > 1, its possibly connected sites, and 
all possible bonds among those sites. The m = 4 case is 
shown in Fig. El To calculate the average clustering co- 
efficient C' m of such a site, we must consider the various 
configurations of long-range bonds among the possibly 
connected sites. The 2 m — 2 potential long-range bonds 
emanating from the original site we divide into two cate- 
gories: the 2 m_1 "shortest" ones (to the sites marked as 
squares in Fig. I17|) . and the remaining 2 m_1 — 2 bonds 
(to the sites marked as triangles in Fig. 1171) . The prob- 
ability for r bonds of the first category and r' bonds of 
the second category is 

/om— 1\ /om— 1 o\ 

?r,r>=[\ )( 2 r , y^'d-pr- 2 — '. 

(Al) 

For a given r and r' , the site has k r y = 2 m + r + r' 
connected sites, so its average clustering coefficient is 

2 m-l 2m -l_ 2 

r=0 r /_ rv r>r i \n,r t r> S-)j^ 

where B r r > is the average number of bonds which actu- 
ally exist among the k r ^ r i sites connected to the original 
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FIG. 17: The open circle on top represents a site with k nn = 
2 m , kid = 2 m — 2, for m = 4, and the figure shows all possible 
neighbors of this site, together with the possible bonds among 
those neighbors. The 2 m sites connected to the top site by 
nearest-neighbor bonds are drawn as filled circles, the 2 m_1 
sites potentially connected to the top site by the shortest long- 
range bonds are drawn as squares, and the other 2 m ~ 1 — 2 
potential neighbors as triangles. 



site. Each of the r bonds of the first category contributes 
two to B r y , as can be seen in Fig. El where there are 
nearest-neighbor bonds connecting every square site to 
two of the 2™ filled circle sites. There are ( r ~t, r ) ways 
of choosing pairs among the r + r' neighbors connected 
to the main site by long-range bonds, but of these pairs, 
only a fraction (2 m — 3)/( 2 2 ~ 2 ) corresponds to possible 
long-range bonds between those neighbors, and of these 
possible bonds on average only a fraction p will actually 
exist. So the total expression for B r ^ r i is 



B r y = 2r + p 



r + r'\ 2 m - 3 



r- 2 ) 



(A3) 



Putting together Eqs. I|A1|) - I|A3() yields the expression for 
C m in Eq. gj). 



2. Average shortest-path length £ n 
for p = and p = 1 

Let us denote the set of sites making up the lattice after 
n construction steps as L n . Then the average shortest- 
path length for L n is defined to be: 



f — 



fin 



N n (N n - l)/2 ' 



where 



fin — ^ dij , 



(A4) 



(A5) 



and dij is the length of the shortest path between sites 
i and j. For the cases p = and p = 1, the lattice 
has a self-similar structure that allows one to calculate 
i n analytically. As shown in Fig. ^] the lattice L n+ i in 
these cases is composed of four copies of L n connected 



(a) 



a 



1, 



at the edges, which we label L„ 
write the sum over all shortest paths fi^+i as 



, 4. We can 



fin 



(A6) 



where A„ is the sum over all shortest paths whose end- 
points arc not in the same L n branch. The solution of 
Eq. (SB is 



S n = 4"- 1 5 1 



n-1 



4™- 



(A7) 



The paths that contribute to A n must all go through at 
least one of the four edge sites (A, B, C, D) at which 
the different L n branches are connected. The analytical 
expression for A„, which we call the crossing paths, are 
found below for p — and p = 1. 



Crossing paths A n for p = 



Let A"'^ denote the sum of all shortest paths with 
endpoints in and L„ . If and meet at 
an edge site, A"-' 3 excludes paths where either endpoint 

is that shared edge site. If and do not meet, 
A" ,/3 excludes paths where either endpoint is any edge 
site. Then the total sum A„ is given by 

A„ = A*- 2 + A 2 ' 3 + A 3 ' 4 + A^ 1 + A^ 3 + A 2,4 

- 2 • 2" +1 . (A8) 

The last term at the end compensates for the overcount- 
ing of certain paths: the shortest path between A and 
C, with length 2™ +1 , is included in both A 2 ' 3 and A 4 ' 1 . 
Similarly the shortest path between B and D, also with 
length 2"+\ is included in both A^ 2 and A 3 ' 4 . 

By symmetry, A,\< 2 = A 2 ' 3 = A 3 ' 4 = A 4 ' 1 and A 4 ' 3 = 
A 2 ' 4 , so that 



A„ = 4A 4 - 2 + 2A 4 - 3 -2-2 



(A9) 
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solution of f n by induction: 

2 „_i 



/„ =J2(2k- 1)2" + J2 2fei/n-i(A) 



fc=i 

2 3n ~ 2 + 2/„_! 



fe=i 



^ 2 fc 2 3(n-fe)-2 + 2 n ~ 1 f 1 



k=0 



2™(2 + 4"), 



(A13) 



where we have used /i = ^i(1) + 2j/ 1 (2) = 4. Substituting 
Eq. jAl~3|) and N n = |(2 + 4") into Eq. jA"lO)) . 



A*' 2 = |2 1+T1 (l + 2 1+2n )(2 + 4") 



(A14) 



FIG. 18: For p — or p — 1, the lattice after n+1 construction 
steps, L n +i , is composed of four copies of L n connected to one 
another as above. The p = 1 case is shown above; for p — 
the horizontal long-range bond is absent. 



A^' 2 is given by the sum 



Proceeding similarly, 
A,\< 3 = £ 4 



A^' 2 = E ^ 



i.i^A 



i#A,D, j=£B,C 

E (diA + 2 n + d, B ) 

+ E ( diD + 2 " + 



E (<^A + <£aj) 



+ 



E 



(A15) 



i.i^A 



(JV„-1) E dw + (JVn-l) E dA > 



2(N n -l) E d iA 



(A10) 



where we have used X^g^t 1 ) ^iA = ~Y^j eL m ^Aj- To find 
y^j C£ W diA, we examine the structure of the hierarchical 

lattice at the nth level. L„ contains v n {m) points with 
diA = m, where 1 < to < 2™, and v n (m) can be written 
recursively as follows: 



The first and second terms are equal and denoted by <?„, 
and the third term is denoted by h n , so that A*' 3 = 
2g n + h n . The quantity g n is evaluated as follows: 

2 «_2 2 n -l-m 

g n = E E I/ n("i)i/„(m')("i + 2 n + m') 

m=l m'=l 
2™ _1 -2 2"~ 1 -l-fc 



f n (m) 



2™ 



if to is odd , 



v n —i(m/2) if to is even . 



(All) 



Expressing X^ e L (1) diA in terms of v n (m), 

2" 

fn = E dlA = E mi/ n(m) . (A12) 



(i) m=l 



Eqs. (|Alip and (|A12|) relate /„ and /„_i, allowing the 



= E E ^-i(fc)^n-i(fc')(2fc + 2" + 2fc') 
fc=i fe'=i 

+ E E ^n-i(A:)2"(2fc + 2" + 2fc'- 1) 
fc=i fe'=i 

2 „_l_ 1 2 T 1 -l_ A , 

+ E E 2 I V„_i(fc')(2/c-l + 2 n + 2fc') 
fc=i fe'=i 

+ E E 2 2n (2fc-l + 2" + 2fc'-l). 
fe=i fe'=i 

(A16) 

The fourth term can be summed directly, yielding 

8 «-i(2«-2)(5-2"-2)/3. (A17) 
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The second and third terms in Eq. (|A16|> are equal and 
can be simplified by first summing over k' , yielding 



2 n - 2 E ^-i(fc)(3-4"-2 Il+2 fc-4£; 2 ). (A18) 



fe=i 



For use in Eq. EED ) T,T=l * v n-i(k) = N n -x - 2, and 
using Eq. (|A13(I . 



2 «-i_ 1 



2 n_ 



fc=l 



fe=l 

■in— 1 //in— 1 



= 2 "- i (4"- i - l)/3. 
Analogously to Eq. ilA13(l , we find 



(A19) 



V k 2 u n ^(k) = l2 2n - 3 {U + 4" - 3n - 3) - 2 r ' 
— ' 9 



fc=i 



(A20) 
(A21) 



With the latter results, Eq. (| A 18(1 becomes 

8 »-i(„23 + 5-4" + 3n)/9. 

With Eqs. (|A17|) and [jA"2l|> . Eq. [jA"l6|) becomes 

g n = 2g n -i + 8™" 1 (-34 - 9 • 2™+ 2 + 25 • 4" + 6n)/9 . 

(A22) 

Using gi = 0, Eq. (|A22|) is solved inductively: 

g n = 2™(164 - 126 • 4™ - 108 • 8™ + 70 • 16™ 

+ 21 -4 n n)/189. (A23) 

All that is left to find an expression for A*' 3 is to evaluate 

2'*-l 

K -2" +1 v n {m)v n {2 n - m) 

m— 1 
2 n -l 



= 2«+ 1 ^, 2 (m) 



m— 1 
2 n- 



= 2 n+1 



^4«+ e wLxw 
fe=i fe=i 

= 16" + 2/i n _i , (A24) 

where the symmetry v n {m) = v n {2 n — m) was used. Us- 
ing hi — 16, Eq. I|A24|I is solved inductively: 

h n = 2™ +3 (8™ - I) /7. (A25) 

From Eqs. (|Al3l) and l(X25|l . 

A i,a = 2 «+! (g _ 18 • 4" + 10 • 16" + 3 • 4 n n) /27 . 

(A26) 

Substituting Eqs. l|AT4|) and lfX26t> into Eq. (JX9J), we ob- 
tain the final expression for the crossing paths A„ when 
p = 0: 

A„ = 2 2+ "(-7 + 3(4 + n)A n + 22 • 16 n )/27 . (A27) 



b. Crossing paths A n for p = 1 
In the p = 1 case, A n is 

A*- 2 + A 2,3 + A 3 ' 4 + A 4,1 + A r \- 3 + A 2,4 - 3 . (A28) 

The last term compensates for the overcounting of the 
shortest path between A and C, with length 2, and the 
shortest path between B and D, with length 1. 

Again by symmetry, A*' 2 = A 3 ' 4 , A 2 ' 3 = A 4 ' 1 , and 
A,\' 3 = A 2 ' 4 , so that 

A„ = 2A, 1 ; 2 + 2A 2 ' 3 + 2A, 1 ; 3 - 3 . (A29) 

We define 

d « 0t = E d iA , 

<c ar ^ E **a, K cai ^ E 



iiA<dio 



diA<d in 



= E < ar ^ E x > ( A3 °) 

diA>dio diA>d iD 

so that c4 ot = < Gar +^ lid + dff r and JV n = A; lcar + N™ id + 
N^. By symmetry 7V£ oar = iV**. Thus, 

A n 2 = £ *i= E (<*iA+<*Aj) 



i 6 i«, je L< 2 > 



i,j^A i.ji^A, d iA <d iD 

E {d iA +d Aj ) 

diA> d iD , djA<dji 

E (d iD + l + d Bj ) 

d ;A><*sD. d jA> d jB 
= £ [(JV„ - ljrfu + < 0t ] 

iei^ 1 ', i#A 

+ E [ W" + ^ id - ^ + <c ar + d « id ] 

iSii, 1 ', i^A 

+ E re car (*D + i)+< car ] 

iSiJ, 1 ', i^A 

: (JV„ ~ l)« Car + dn" 1 ) + {K Bal + - 

+ (K lcar + K nid - i)4? r + Ar^Cd" 631 + < lid ) 



+ JV£ ear ((C ar + K lcar ) + K lcar < car ■ 



(A31) 



The horizontal long-range bond does not affect A 2 ' 3 , so 
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that Eq. (|JTp)l still holds, A 2 - 3 = 2(N n - l)c^ ot . Finally, 
a^ 3 = E d a = E (d lD + i + d Cj ) 

i^tA,D, j^B,C i^A,D, j^B,C, d jB >d jC 

+ E {dio + l + dBj) 

iGii 1 ', J6ii 1 3) , i^A,D 
j^B,C, d iA >d iD , dj B <dj C 

+ E {d iA +l + d Bj ) 

j'6-ti 3 ', i^A,D 
j^B,C, d iA <d iD , d jB <d jC 

E [(K lcar -i)(rf lD + i) + < car ] 

+ E [(wr+wf-iife + i) 

diA>diD 

+ E [(ivr ar +K nid -i)(^A + i) 







= (A^ Car -l)(<4 0t -l + ^n- 


2) + (N n - 2)< car 


+ (N™ al + N™ id - l)(<C ar 


+ C M + W£ eal + M - 1) 


+ (A r ; lcar + N™ id - 


i)« car + C id ) 


+ (7V£ car + N™ id - l)« car 


+ ^ car - 1) 


+ (Af; lcar - l)« car 


+ C id )- (A32) 


Having A*' 2 , A 2 ' 3 , and A*> 3 in 


terms of the quantities in 



Eq. (|A30|1 . the next step is to explicitly determine these 
quantities. 

We consider a site i € and the shortest-path dis- 
tances to the edges, diA and diD- If the site was added 
to the lattice at the nth construction step, the values 
of diA and dio do not change at subsequent steps, since 
the shortest path to the edge sites is always along the 
bonds added earliest. We see this in Fig. where the 
sites are labeled by the ordered pairs diD, diA, for the 
first three construction steps. We denote by m , the 
number of sites added at the nth construction step which 
have diA = m, dn> = m! . Since A and D are connected 
by a long-range bond, m' and m can differ by at most 
1. Thus for a given m there are three categories of sites 
added at step n, respectively numbering m+i> a m ra> 
and a" +lm . By symmetry a™ +1 m = a™ m+1 . The m, 
mf values of sites added at step n depend on the neigh- 
boring sites, which were added at previous construction 
steps. For example, there are 2 ra sites added at the nth 
step (n > 2) which are nearest- neighbors of site A, so 
these new sites have m = 1, m' = 2, giving a™ 2 = 2". 
Sites with m = 1, m' = 2 will in turn get neighbors with 
m = 2, m' = 3 in subsequent steps. The relationship 
between a 2 3 and a\ 2 for k < n is 

n-2 

Oa,3 = E 2" +1 ^ fc af, 2 = 2" +1 (n - 3) . ( A33) 

fe=2 



1,1 




FIG. 19: The first three construction steps of the lattice with 
p = 1, with the sites labeled by ordered pairs denoting the 
shortest-path distance to the left and right edge sites. 



Similarly, 

<4 = E 2 n+1 ~ k <3 = 2" +1 (« - 4)(n - 5) . (A34) 

Since sites with distances m, m + 1 do not appear before 
the construction step n = 2m, the sum over a 2 3 starts 
at k = 4. Proceeding in this manner, for general m > 1 
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and n > 2m, 



n-2 

E2«+i-fc fc 
m — l,r. 

fe=2(m-l) 

2 m - 1 2"(n-m- 1)! 



(m- l)!(n-2m)! 



(A35) 



The value of Og i is 1 for n = and for n > 0. Analo- 
gously, for general m > 2 and n > 2m — 1, 



2 m - 1 2 n (n-m-l)l 
(m-2)!(ra-2m+l)! 



(A36) 



2 r 2 „+3 _ 12 . 2 n„ + ( 3n + 7 ) 4 n + jl 



2_ + y . 4 n + 3n ( 2 2+n + 4 ") + ]_] 



(A37) 



where [^J denotes the largest integer < a; and the dif- 
ferent results for n even and odd are given consecutively, 
and 



incar i i\rncar 



n L«'/2J 

< ar = E E + 

n'=l fc=l 



Substituting the results of Eq. lA37jl into Eqs. I|A31I) - 
31J|,forp = l, 



The value of a" j is 2 for n = 1 and for n > 1 . 
Thus we obtain the quantities in Eq. (|A30|) . 

n |n'/2J 
jyncar _ 1m,m + l = { ^ 



1+ §(2 n +4 n - 2) 

^ + §(2" - 2)(2" + 1) 



"' 1 |(2 n + l) 2 



„ L(«'+i)/2J [2,™ n a 

N™ d = E E < J " 

n' = l fc=l 
n L«'/2J 

< car = E E 

n'=l fc=l 
n L(«' + 1)/2J 

■t'r = E E 



= J^( 2n + 2 )( 2n +3-2"n-l) 
1 A(2 n - 2)(2 n +3 ■ 2 n n + 1) 



k 



A„ = — [-23 - 8(-2) n + 8 • 4" + (104 + 48n)16 n ] . 

(A39) 

Substituting Eqs. (|A27j) and (|A39fl for A m into 
Eq. iJXTJl, and using Si = 8, 7 for p = 0, 1, 



_ f [98 + 27 • 2™ + (42 + 21n)4" + 22 ■ 16™] 
" _ j i [23 + 4(-2)" + (44 + 6n)4™ + (10 + 12n)16 n ] 

(A40) 

Eq. I|A40|) in Eq. (|A4|) yields the analytical expressions 
for £ n in Eqs. and ©. 
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